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>■ , Given a completely integrable system, we associate to any connection on a fiber bundle 
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in invariant tori over a parameter manifold the classical and quantum holonomy operator 
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(generalized Berry's phase factor), without any adiabatic approximation. 



O : I. INTRODUCTION 

At present, holonomy operators in quantum systems attract special attention in connec- 
tion with quantum computation (see, e.g., Refs. [1-3]). They exemplify the non-Abelian 
generalization of Berry's geometric phase by means of driving a finite level degenerate eigen- 
state of a Hamiltonian over a parameter manifold. The key point is that a geometric phase 
depends only on the geometry of a path executed and, therefore, provides a possibility to 
perform quantum gate operations in an intrinsically fault-tolerant way. The problem lies 
in separation of a geometric phase factor from the total evolution operator without us- 
ing an adiabatic assumption. Firstly, holonomy quantum computation implies exact cyclic 
evolution, but exact adiabatic cyclic evolution almost never exists. Secondly, an adiabatic 
condition requires that the evolution time must be long enough. 

A nonadiabatic Abelian phase was discovered by Aharonov and Anandan who con- 
sidered a loop in a projective Hilbert space instead of a parameter space. 4 Non-Abelian 
generalization of the Aharonov-Anandan phase has been studied under rather particular 
assumption. 5 . Moreover, a non-Abelian Aharonov-Anandan phase fail to be separated from 
the dynamic one in general. Recently, several schemes using the Aharonov-Anandan phase 
were proposed for nonadiabatic geometric gates. 6-8 
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In a general setting, let us consider a linear (not necessarily finite-dimensional) dynam- 
ical system d t ip = Sip whose linear (time-dependent) dynamic operator S falls into the 
sum 

S = S + A = S + A a d t C, (1) 

where is a function of time taking its values in a finite-dimensional smooth real param- 
eter manifold E coordinated by (<J a ). Let assume that: (i) the operators S (t) and A(t') 
commute for all instants t and t', and (ii) the operator A depends on time only through 
Then the evolution operator U(t) can be represented by the product of time-ordered 
exponentials 



U(t) = U (t) o U^t) =Texp 



where the second one is brought into the ordered exponential 
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U^t) = Texp 



j A a {i(t'))d t e(t')dt' 
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= Texp 



A a (a)da a 



(3) 



along the curve £[0,t] in the parameter manifold E. It is a nonadiabatic geometric factor 
depending only on the trajectory of the parameter function £. Accordingly, A is a holonomy 
operator. The geometric factor (3) is well defined if A a da a is an Ehresmann connection on 
a fiber bundle over a parameter manifold E. Then this factor is a displacement operator 
along an arbitrary curve £[0, t] C E. 

A problem is that the above mentioned commutativity condition (i) is very restrictive. 
Moreover, it need not be preserved under time-dependent transformations. 

For instance, let us consider a Hamiltonian system of dynamic variables (q,p). Written 
with respect the initial data coordinates (qo,po), its Hamiltonian 7Y(g ,R)) vanishes. Given 
these coordinates (qo,Po), let one can introduce a perturbed Hamiltonian H.£(qo,Po,£(t)) 
which depends on parameter functions £(t) and generates a holonomy operator A (1). Then 
the evolution operator of the perturbed Hamiltonian system reduces to the geometric factor 
(3). Relative to the original variables (q,p), a Hamiltonian of this perturbed Hamiltonian 
system is 

W = H(q,p,t) + Hz(q (t,q,p),p(t,q,p),t(t)). 

However, the corresponding evolution operator does not fall into the product (2) because a 
Hamiltonian H is not a function under time-dependent transformations and, consequently, 
the Poisson bracket {H, H{\ with respect to original variables (q,p) need not vanish. 

Nevertheless, basing on this example, we can essentially extend the class of dynamical 
systems admiting a nonadiabatic geometric phase. We aim to describe dynamical systems 
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where the commutativity condition (i) is not satisfied, but a part of dynamic variables is 
driven only by a holonomy operator. These are completely integrable Hamiltonian systems. 

Let us consider a completely integrable Hamiltonian system (henceforth CIS) of m 
degrees of freedom around its invariant tori T m . We show that, being constant under an 
internal evolution, its action variables are driven only by a perturbation holonomy operator 
A which can be associated to an arbitrary connection on a fiber bundle 

ExT^S. (4) 

This holonomy operator is defined with respect to the initial data action-angle coordinates 
without any adiabatic approximation. Then we return to the original action-angle coordi- 
nates. The key point is that both classical evolution of action variables and mean values 
of quantum action operators relative to original action-angle coordinates are determined in 
full by the dynamics of initial data action and angle variables. 

The plan of the paper is as follows. Section II addresses classical time-dependent CIS. 
The key point is that any time-dependent CIS of m degrees of freedom is extended to an 
autonomous CIS of m + 1 degrees of freedom 9-11 and, as a consequence, can be provided 
with action-angle variables around a regular instantly compact invariant manifold. 10 ' 11 

In Section III, we introduce the holonomy operator in a classical CIS by use of the 
fact that a generic Hamiltonian of a mechanical system with time-dependent parameters 
contains a term which is linear both in momenta and the temporal derivative of a parameter 
function. 12 ' 13 This term comes from a connection on the configuration space of the system 
fibered over a parameter manifold. 

Section IV is devoted to geometric quantization of a CIS with respect to the angle 
polarization. This polarization leads to the Schrodinger representation of action variables 
in the separable Hilbert space of smooth complex functions on y m . 10 ' 14 We show that this 
quantization both with respect to the original action-angle variables and the initial data 
action-angle variables is equivalent. 

In Section V, the classical holonomy operator of Section III is quantized with respect 
to the initial data action-angle variables. 

The symbols J and |_ below stand for the left and right interior products of multivector 
fields and exterior forms, respectively. 

Let us recall that, given a fiber bundle Y — > X coordinated by (x x ) y l ) ) a connection K 
on Y — > X is defined by a tangent-valued form 

K = dx x ® (<9 A + k\di) 

on Y. 15 A connection on a fiber bundle Y — * X is said to be an Ehresmann connection if, 
given an arbitrary smooth curve £([0, 1]) C X, there exists its horizontal lift through any 
point of Y over £(0). 
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Let X be a real axis R provided with the Cartesian coordinate t possessing transition 
functions t' = t+const. A connection K on a fiber bundle Y — > R is uniquely represented by 
a vector field K on Y such that K\dt = I. 12 This is the case of time-dependent mechanics. 

II. CLASSICAL TIME-DEPENDENT CIS 

Recall that the configuration space of time- dependent mechanics is a fiber bundle Q — > R 
over the time axis R. Let it be equipped with the bundle coordinates (t, q k ), k = 1, . . . , m. 
The corresponding phase space is the vertical cotangent bundle V*Q of Q — > R endowed 
with the induced coordinates {t,q k ,Pk) relative to the holonomic coframes {dq k }. 12,16 The 
cotangent bundle T*Q of Q — > R plays a role of the homogeneous phase space of time- 
dependent mechanics. It is equipped with the induced coordinates (t, q k ,p,Pk) relative to the 
holonomic coframes {dt,dq k }. With respect to this coordinates, the canonical symplectic 
form and the corresponding Poisson bracket on T*Q read 

Q = dp Adt + dp k A dq k , 

{/, /'}r = 9 P M/' - 9 t /9 p /' + 9 fc /5 fc f - 4/9 fc r, /, /' e C°°(T*Q). 

There is the one-dimensional trivial afline bundle 

C : T*Q -> \/*g. (5) 

As a consequence, the phase space V*Q of time-dependent mechanics is provided with the 
canonical Poisson structure 

{/, f'}v = d k fd k f - d k fd k f, /, /' G C°°(V*Q), (6) 

given by the relations 

C{f, f'}v = {Cf, Cfh, f, f e C°°(V*Q). 

The corresponding Poisson bivector on V*Q reads wy = d k A d k . 

A Hamiltonian of time-dependent mechanics is defined as a global section 

h:V*Q^T*Q, poh = -H(t,q j , Pj ), (7) 

of the afline bundle ( (5). 12,16 Given the pull-back form h*Q, the relations 7_ffJ<it = 1, 
•jH\h*Q = define a unique Hamilton vector field 

lH = d t + d k Hd k - d k Hd k (8) 
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on V*Q and the corresponding Hamilton equations 

q k = d k H, p k = -d k H. (9) 

Note that, given a connection T = dt + T\di on Q — » R, any Hamiltonian 7i (7) admits 
the decomposition 7i = PiT\ + H where TL is a function on V^Q- 

An integral of motion of the Hamilton equations (9) is a smooth real function F on V*Q 
whose Lie derivative 

L, H F = lH \dF = d t F + {H,F} v 

along the Hamilton vector field 7# (8) vanishes. A time-dependent Hamiltonian system 
of m degrees of freedom is a CIS if there exist m independent integrals of motion {F k } in 
involution with respect to the Poisson bracket {, } v (6). Their Hamiltonian vector fields 

tf. = - Wv [dFi = d k F t d k - d k Fd k 

and the Hamilton vector field 7# (8) generate a smooth regular distribution on the phase 
space V*Q and the corresponding foliation of V*Q in invariant manifolds. 

One can associate to any time- dependent CIS on V*Q an autonomous CIS on the 
homogeneous phase space T*Q as follows. 

Given a Hamiltonian h (7), let us consider an autonomous Hamiltonian system on the 
symplectic manifold (T*Q, Q) with the Hamiltonian 

H = d t \(~-Ch*~)=p + H. 

Its Hamiltonian vector field 

7r = dt - d t Hd p + d k Hd k - d k Hd k (10) 
is projected onto the Hamilton vector field 7^ (8) on V*Q so that 

C(W) = {h,C7}t, feC°°(v*Q). 

An immediate consequence of this relation is the following. 

(i) Given a time-dependent CIS (if; Fk) on V*Q, the Hamiltonian system {H, (*F k } on 
T*Q is a CIS. 

(ii) If M C V*Q is an invariant manifold of the time-dependent CIS {H;F k }, then 
h(M) C T*Q is an invariant manifold of the homogeneous CIS (H, (*F k ). 

Hereafter, let the Hamilton vector field 7^ (8) be complete, i.e., the Hamilton equations 
(9) admit a unique global solution (a trajectory of 7^) through every point of the phase 
space V*Q. The trajectories of 7^ define a trivial bundle V*Q — > VqQ over the fiber VqQ 
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of V*Q — > R at t = 0. Then any invariant manifold M of {H; F^} is also a trivial bundle 
M = lx M over M = M n VqQ. 

If M is compact, one can introduce action- angle coordinates around an invariant ma- 
nifold M by use of the action-angle coordinates around the invariant manifold h(M) of 
the corresponding autonomous CIS on T*Q. 10 Namely, h(M) has an open neighbourhood 
which is a trivial bundle 

U' = V' x R x T m -> V x R -> V' (11) 

over a domain V' C R m+1 with respect to the action-angle coordinates (Iq, Ii,t,4> 1 ). Here- 
with, the following holds, (i) Iq = H. (ii) The integrals of motion (*Fk depend only on the 
action coordinates Jj. (iii) The symplectic form Q on U' reads 

Q = dl A dt + dli A dcji . 

The symplectic annulus U' (11) inherits the fibration structure (5) over the toroidal domain 

U = V x R x T m , V cR m . (12) 

Coordinated by (Jj,t,0*) and provided with the Poisson structure (6), the toroidal domain 
(12) is a phase space of the time-dependent CIS (H; Fj) around its instantly compact 
invariant manifold M. Since H = J , the Hamilton vector field (10) is 7t = <9t, and so is its 
projection 7^ (8) onto U. Hence, the above-mentioned action-angle coordinates (h,t, 
are the initial data coordinates. 

These action-angle coordinates are by no means unique. Let 7i be an arbitrary smooth 
function on R m . Then the canonical transformation 

I'o^Io-Hilj), Il = h, t' = t, i P i = f + td l H{I 3 ) (13) 

gives new action-angle coordinates corresponding to a different trivialization of U' (11) (and 
U (12)). Accordingly, the Hamilton vector field 7# takes the form (8), and the Hamilton 
equations (9) read 

v? fc = o k n, i k = 0. 

These are the Hamilton equations of an autonomous CIS with a time-independent Hamil- 
tonian H on the toroidal domain U (12). 

III. CLASSICAL HOLONOMY OPERATORS 

The phase space of a Hamiltonian system with time-dependent parameters is a compos- 
ite fiber bundle n — > E x R — > R, where n — > E x R is a symplectic bundle and Sxl^R 
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is a parameter bundle whose sections are parameter functions. 12,13,17,18 In the case under 
consideration, all bundles are trivial and their trivializations hold fixed. Namely, the phase 
space is the product 

n = Sx(/ = Sx(^xlxr)-4SxK^R, 

equipped with the coordinates (a a , I k , t, (f> h ). Let us suppose for a time that parameters are 
also dynamic variables. The phase space of this system is the fiber bundle 

n' = t*e x[/^sxixr 

coordinated by (cr a , a a , I k , t, <p k ). A generic Hamiltonian of such a system is 

Hx = <r a Z? + I k (A k + A*E«) + H(aP, I 3 , t, <fP), (14) 

where 

d t + + (A k + A^)d k 

is a composite connection on the fiber bundle £ x R x T rn — > R generated by a connection 
d t + S"(9 a on the parameter bundle Sxl->K and a connection 

A = dt <g> (d t + A k d k ) + da a <g) (0 a + A%) (15) 

onExlx T m — > £ x R. 12,13,18 Then a Hamiltonian system with a fixed parameter function 
a a _ i s characterized by the Hamiltonian 

= i k \A k t (t, <fp) + a*(^, t, <mr] + ^, * , <p j ) (is) 

on the pull-back bundle U = £*n (12). 

Let (Ik,t,(f) k ) be the initial data action-angle coordinates of a time- dependent CIS. Its 
Hamiltonian 7i with respect to these coordinates vanishes. Therefore, we can introduce a 
desired holonomy operator by the appropriate choice of the connection A (15). Let us put 
A k = and assume that coefficients A k a are independent of time, i.e., the part 

A s = da a ® (d a + A k a d k ) (17) 

of the connection A (15) is a connection on the fiber bundle (4). Then the Hamiltonian of 
a perturbed CIS reads 

n^ = i k A k a (^,^)dtC. (is) 

Its Hamilton vector field (8) is 

lH = dt + KdtCdi - IAAldtCd 1 - (19) 
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It leads to the Hamilton equations 

^ = A^(t),0^r, (20) 

d t I t = -I k d t A k a mA l )dtC- (21) 

Note that 

V*A S = da a <g> (d a + A^di - hdiAffi (22) 

is the lift of the connection (17) onto the fiber bundle £ x (V x T m ) — > £, seen as a 
subbundle of the vertical cotangent bundle V*(Y> x T m ) = £ x T*T m of the fiber bundle 
(4). It follows that any solution l (t) of the Hamilton equations (20) - (21) (i.e., an 
integral curve of the Hamilton vector field (19)) is a horizontal lift of the curve £(£) C E 
with respect to the connection V*A S (22). i.e., hit) = Ii(£(t)), ^(t) = </>(£(t)). Thus, the 
right-hand side of the Hamilton equations (20) - (21) is the holonomy operator 

A = (A^C , -hdiA k a d t C) (23) 

(cf. (1) where §0 = 0). It is not a linear operator, but the substitution of a solution 0(£(t)) 
of the equation (20) into the Hamilton equation (21) results in a linear holonomy operator 
on the action variables Jj. 

Let us show that the holonomy operator (23) is well defined. Since any vector field d 
on R x T m such that d\dt = 1 is complete, the Hamilton equation (20) has solutions for 
any parameter function It follows that any connection A^ (17) on the fiber bundle (4) 
is an Ehresmann connection, and so is its lift (22). Therefore, any curve £([0, 1]) C E can 
play the role of the parameter function in the holonomy operator A (23). 

Now, let us return to the original action-angle coordinates (Ik,t,(f k ) by means of the 
canonical transformation (13). Relative to these coordinates, the perturbed Hamiltonian 
reads 

W = i k A k Mt), & - tom^MCit) + nil,), 

and the Hamilton equations (20) - (21) take the form 

W = amy,) + Aim, v l - td'n^mcit) 

-thdwn^dsA^t), <p l - td l n{i 3 ))d t C{t), 
d t ii = -i k d t A k a m, v l - td l n(i 3 ))d t c(t). 

Their solution is I^(t)), ^(t) = 4 (£(t)) + t<9^(/,(£(t))) where I^(t)), 0*(£(t)) is a 
solution of the Hamilton equations (20) - (21). It is readily observed that the action 
variables l k are driven only by the holonomy operator, while the angle variables ip l have a 
nongeometric summand. 
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Let us emphasize that, in the construction of the holonomy operator (23), we did not 
impose any restriction on the connection (17). Therefore, any connection on the fiber 
bundle (4) generates a holonomy operator in a CIS. However, a glance at the expression 
(23) shows that this operator becomes zero on action variables if all coefficients A k of the 
connection As (17) are constant, i.e., As is a principal connection on the fiber bundle (4) 
seen as a principal bundle with the structure group T m . 

IV. QUANTUM CIS 

There are different approaches to quantization of CISs. 19 Their geometric quantization 
was studied at first with respect to the polarization spanned by Hamiltonian vector fields 
of integrals of motion. 20 For example, the well-known Simms quantization of the harmonic 
oscillator is of this type. In this approach, the problem is that the associated quantum al- 
gebra includes afline functions of angle coordinates which are ill defined. As a consequence, 
elements of the carrier space of this quantization fail to be smooth, but are tempered 
distributions. In recent works, 10 ' 14 we have developed a different variant of geometric quan- 
tization of CISs by use of the angle polarization spanned by almost-Hamiltonian vector 
fields d k of angle variables. This quantization is equivalent to geometric quantization of the 
cotangent bundle T*T m of a torus T m with respect to the vertical polarization. The result 
is as follows. 

Given an autonomous CIS on a symplectic annulus 

p = v xT m , n P = dlt A dip* 

equipped with the action-angle coordinates (ij, ip 1 ), its quantum algebra A with respect to 
the above mentioned angle polarization consists of afline functions 

/ = aV')4 + % J ) 

of action coordinates I k . They are represented by self-adjoint unbounded operators 

/ = - ta k d k - l -8 k a k - a k \ k + b (24) 

in the separable pre-Hilbert space of complex half-forms on T m . If coordinate transforma- 
tions of T m are only translations, this space can be identified with the pre-Hilbert space 
C 00 (T m ) of smooth complex functions on T m . Different tuples of real numbers (Ai, . . . , A m ) 
and (X[, . . . , X' m ) specify inequivalent representations (24), unless \ k — \' k e Z for all 
k — 1, . . . , m. These numbers come from the de Rham cohomology group iJ 1 (T m ) = M m . 

In particular, the action operators (24) read I k = —id k — \ k . They are bounded. By 
virtue of the multidimensional Fourier theorem, an orthonormal basis for C°°(T m ) consists 
of functions 

VV)^) = exp[m r y? r ], (n r ) = (m, . . . , n m ) e Z m . (25) 
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With respect to this basis, the action operators are brought into countable diagonal matrices 



= (nk - A fe )VV)> ( 26 ) 

while functions a k ({p) are decomposed in Fourier series of the functions ip( nr ), which act on 
C°°(T m ) by the law 

VVoV'K) = V'K+o- (27) 

It should be emphasized that a I k 7^ a k I k 7^ 1/^. 

If a Hamiltonian Ti.(Ik) of an autonomous CIS is an analytic function on R m , it is 
uniquely quantized as a Hermitian element 7i{Ik) = 7~L(Ik) of the enveloping algebra of A. 
It is a bounded self-adjoint operator with the countable spectrum 

H(I k )lp(n r ) = £(n r )VV-)> E (n r ) = ^( n k - Afc), U k E (n r ). (28) 

In order to quantize a time-dependent CIS on the Poisson toroidal domain (U, {, }y) (12) 
equipped with action-angle coordinates (Jj, t, ip l ), one may follow the instantwise geometric 
quantization of time-dependent mechanics. 21 As a result, we can simply replace functions 
on T m with those on M x T m . w Namely, the corresponding quantum algebra A C C°°(U) 
consists of affine functions 

f = a k (t,if^)I k + b(t,^) (29) 
of action coordinates I k represented by the operators (24) in the space 

E = C°°(1R x T m ) (30) 

of smooth complex functions ip(t, ip) on R x T m . This space is provided with the structure 
of the pre-Hilbert C°°(lR)-module with respect to the nondegenerate C°°(IR)-bilinear form 

—J J ^'d m <p, i\),i>' G C°°(R x T m ). 

Its basis consists of the pull-back onto R x T m of the functions VW) (25). 

This quantization of a time-dependent CIS is extended to the associated homogeneous 
CIS on the symplectic annulus (U',Q) (11) by means of the operator I Q = —id t in the 
pre-Hilbert module E (30). Accordingly, the homogeneous Hamiltonian H is quantized as 
H = —id t + 7~L. The corresponding Schrodinger equation is 

H0 = -id t ip + Hip = 0, 4>EE. (31) 

For instance, the quantum Hamiltonian of the original autonomous CIS is 

U = -idt + H{Ij). 
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Its spectrum H/0(„ r ) = E(n r )^< nr ) relative to the basis {ip(n r )} for E (30) coincides with that 
of the autonomous Hamiltonian (28). The Schrodinger equation (31) reads 

= -idtip + H{-id k - A fe )V = 0, ip e E. 

Its solutions are the Fourier series 

i> = J2 B (n r ) exp[-**^(nr)]^(n r )> B (n r ) £ C. 

(n r ) 

Now, let us quantize this CIS with respect to the initial data action-angle coordinates 
(Ii,(f) 1 ). Its quantum algebra Aq C C°°(U) consists of affine functions 

f = a%<fS)I k + b(t,<fS). (32) 

The canonical transformation (13) provides an isomorphism between Poisson algebras A 
and Ao- Functions / (32) are represented by the operators / (24) in the pre-Hilbert module 
Eq of smooth complex functions ^f(t, (j)) on R x T m . Given its basis ^/( nr )(0) = [m r r ], the 
operators Ik and ip(n r ) take the form (26) and (27), respectively. The Hamiltonian of a 
quantum CIS with respect to the initial data variables is H = —id t . Then one easily 
obtains the isometric isomorphism 

i?(Vv)) = exp[^£ K) ]* K) , W) W)> = (W), (33) 

between the pre-Hilbert modules E and E which provides the equivalence 

T i = R- 1 T i R, $ M = R-% nr) R, U = R~ l U R (34) 

of the quantizations of a CIS with respect to the original and initial data action-angle 
variables. 

V. QUANTUM HOLONOMY OPERATORS 

In view of the isomorphism (34), let us first construct a holonomy operator for a quantum 
CIS (Aq, H ) with respect to the initial data action-angle coordinates. Let us consider the 
perturbed homogeneous Hamiltonian 

H ? = H + H x = J + d t C(t)A k Mt), <P)h 

of the classical perturbed system (18). Its perturbation term Hi is of the form (29) and, 
therefore, is quantized by the operator 

Hx = -id t CK = -td t C[A h a d k + \d k (k k a ) - i\ k k k a \. 
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The quantum Hamiltonian = H + Hx defines the Schrodinger equation 

dp + d t C[A% + i^(Aj) - iA fe A*]* = 0. (35) 

If a solution exists, it can be written by means of the evolution operator which reduces to 
the geometric factor U\ (3). The latter can be viewed as a displacement operator along the 
curve £[0, l]cS with respect to the connection 

A s = da a (d a + A a ) (36) 

in the C°°(£)-module C°°(S x T m ) of smooth complex functions on £ x T m . 13 ' 15 ' 18 ' 22 Let 
us study the existence if this displacement operator. 

Given a connection A s (17), let $*(t,0) denote the flow of the complete vector field 
d t + Aj,(£, (f))d t £ a di on R x T m . It is a solution of the Hamilton equation (20) with the 
initial data 0. We need the inverse flow ($~ 1 )*(t,0) which obeys the equation 

Let \l/o be an arbitrary complex half-form \l/o on T m possessing identical transition functions, 
and let the same symbol stand for its pull-back onto R x T m . Given its pull-back 

(* -1 r*o = det (^§P^ ' M*-\t, 0)), (37) 
it is readily observed that 

^ = ($-i)*^ o exp[a fc fc ] 

obeys the Schrodinger equation (35) with the initial data \Ev This function is well defined 
only if all the numbers equal or ±1/2. Note that, if some numbers are equal to 
±1/2, then \l/ exp[i\ k (J) k ] is a half-density on T m whose transition functions equal ±1, i.e., 
it is a section of a nontrivial metalinear bundle over T m 

We thus observe that if equal or ±1/2, then the displacement operator always 
exists and A = iHi is a holonomy operator. A glance at the action law (27) shows that 
this operator is infinite-dimensional. 

For instance, let As (17) be the above mentioned principal connection, i.e., A^ =const. 
Then the Schrodinger equation (35) where X k = takes the form 

W(t, V) + d t C(t)A k a d k ^(t, <jt) = 0. 

Its solution (37) is 

$(^) = * (^-r(t)-noK). 
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The corresponding evolution operator reduces to Berry's phase multiplier 



U^ (nr) = exphm^r (t) - £ a (0))A£)]¥ (Bp) , 



rtj E (n r ). 



It keeps the eigenvectors of the action operators Jj. 

In order to return to the original action-angle variables, one can employ the morphism 
R (33). The corresponding Hamiltonian reads H' = i? _1 H^i?. The key point is that, due 
to the relation (34), the action operators ij have the same mean values 



with respect both to the original and the initial data action-angle variables. Therefore, 
these mean values are defined only by the holonomy operator. 



We have shown that any CIS around its compact invariant manifold admits a perturba- 
tion dependent on parameters by means of holonomy operator associated to a connection 
on the fiber bundle (4). 

Since action variables are driven only by a holonomy operator, one can use this operator 
in order to perform a dynamic transition between classical solutions or quantum states of 
an unperturbed CIS by an appropriate choice of a parameter function £. The key point is 
that this transition can take an arbitrary short time because we are entirely free with time 
parametrization of £ and can choose it quickly changing, in contrast with slowly varying pa- 
rameter functions in adiabatic models. For instance, one can choose £ a step function, then 
its time derivative is a 5-function of time. This fact makes nonadiabatic holonomy opera- 
tors in CISs promising for several applications, including classical and quantum scattering 
in integrable Hamiltonian systems, 23 quantum control operators, 24 ' 25 and the above men- 
tioned quantum computation. It also looks attractive that quantum holonomy operators in 
CISs are essentially infinite-dimensional, whereas both the existent quantum control theory 
and the theory of quantum information and computation 26 involve only finite-dimensional 
operators. 



* = R(ri>), 
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